Stability of the self-phase-locked pump-enhanced singly resonant parametric oscillator 
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Steady-state and dynamics of the self-phase-locked {3lu — > 2uj,ui) subharmonic optical parametric 
oscillator are analyzed in the pump-and-signal resonant configurations, using an approximate ana- 
lytical model and a full propagation model. The upper branch solutions are found always stable, 
regardless of the degree of pump enhancement. The domain of existence of stationary states is found 
to critically depend on the phase-mismatch of the competing second-harmonic process. 
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A new class of subharmonic (frequency divide-by-three, 
or 3 1) optical parametric oscillators (OPOs), namely 
self-phase-locked (SPL-OPOs), has aroused lately much 
interest from both the experimental J^^ J] and theo- 
retical d, ^ 1^ standpoints. Consider a 3 1 OPO 
pumped at an angular frequency ujp = 3uj (generating 
signal and idler waves at tOg = 2a; and tOi = lo) con- 
taining a second nonlinear crystal that is phase-matched 
for the competing second-harmonic (degenerate down- 
conversion) of the idler (signal) waves. The loss-free 
X^^) medium of total length L = Li -\- L2 consists of a 
dual-grating periodically-poled (PP) crystal comprising 
a first section of length Li perfectly phase-matched for 



the Xopo 



DRO/TROs (subject to mode pair instabilities) or SROs 
(requirement of high pump thresholds), it is interesting 
to check how the dynamics of the doubly/triply reso- 
nant SPL-OPOs is affected when only one of the sub- 
harmonic waves oscillates with a varying pump enhance- 
ment, and to compare with the behaviour of the purely 
idler-resonant (SPL-IRO) case treated by Longhi in the 
mean-field approximation and neglecting pump deple- 
tion 0. 

The analysis starts with the solutions of the reduced 
propagation equations for the normalized slowly vary- 
ing field envelopes Aj{t', Z) = giLiNj{t', Z) (j = p, s, i) 
throughout the dual-section medium. The complex am- 



'(-3w;2w,w) down-conversion, followed plitudes are scaled such that their intensities Ij = \Aj 



by a second section of length L2 nearly phase-matched 
for the competing Xshg = x''^'' (~2a;; w, w) process with 
a wavevector mismatch Afc = fc2cj — 2fci^ ^ 0. Due to 
the mutual self-injection of the subharmonic waves, the 
dynamics of the signal- and-idler resonant devices was 
recently shown to depart from that of a conventional 
{L2 = 0) non-degenerate OPO Q. The main differ- 
ence is that SPL-OPOs are characterized by an intensity 
bistability (^b-critical bifurcation) whatever the config- 
uration [1, 0, ■ Secondly, each intensity state of the 
stable branch can take 3 possible deterministic phase val- 
ues equally spaced by 27r/3, while conventional OPOs are 
subject to a stochastic phase diffusion process 0. Exper- 
imentally, 3-;- 1 SPL-OPOs are investigated in frequency 
metrology (Ij, i?,, ^ , Fourier synthesis of attosecond pulse 
0, transverse pattern formation 0,0 and potentially for 
new features in squeezed states of light. Finally, owing to 
the much lower pump intensity requirement, SPL-OPOs 
should be also suited for the first experimental evidence 
of a Hopf bifurcation in cw OPOs, predicted to only occur 
in triply-resonant OPOs under extreme detuning condi- 
tions 

In this brief report, I provide an analysis of the pump- 
enhanced singly resonant device (SPL-PRSRO) which 
was only over-viewed in the conclusion of Ref. 3] . Bear- 
ing in mind that PRSROs are easier to implement than 



are proportional to the number of photons \Nj\'^ m j- 
th mode, times the square of the small-signal parametric 
gain giLi (<C 1) of the 3-^1 process {gi oc Xopo) I^- In 
the phase-retarded time frame {t = t' — nZ /c; z = Z/ Li 
or z' = Z/L2), where n is the average index of refraction, 
the 3 plane-wave equations for < z < Li/Li = 1 are 



dzAai — iA„A, 



(1) 



with the initial condition Ai(z = 0) = (non-resonant 

< L2/L2 = 1), the 



idler). In the SHG section (0 < z 
subharmonic amplitudes evolve as 

iSA^e^p (+i2Cz'), 
iSA,A* exp {~i2^z') 



dz'As = 



(2) 
(3) 

.(^ = 1), 
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with the initial conditions As^i{z' = 0) = As^ty 
while the pump amplitude keeps its value at z ~ 1 
throughout the second section {dz>Ap — 0). The pa- 
rameter ^ = AfcL2/2 is the phase mismatch of the com- 
peting SHG process. The nonlinear coupling parame- 
ter S = g2L2/giLi is the ratio of the SHG to OPO 

(2) 

small signal gains (32 oc Xshg)- ^^s expression reduces 
to 5 ~ {y/Lii/Vi for a 3 ^ 1 OPO employing a PP 
material |3| . The time-dependent cavity dynamics is ob- 
tained from an iterative mapping of the resonating field 
amplitudes at z' = 1 and at a time t to their values at 
z — after one roundtrip time r of the ring-type cavity, 

Aj{t + T,z = 0) = rj cxp{iAj)Aj{t, z' = 1) + A„,(4) 

with j = p,s and with rj being the (real) amplitude re- 
flectivity from z' = 1 back to z = 0. The constant input 



2 



field Ain stands for the driving pump field and is null 
for the signal wave. The rj's are related to the cavity 
loss parameters k^'s by Hj — I — rj. For the resonat- 
ing subharmonic, it will be always assumed that Kg <^ 1, 

~ 1. The amplitude loss parameter Ks is then related 
to the cavity finesse by = tt/ks and to the cavity half 
linewith by 2TT"fs = i^s/t. The phase factors Aj corre- 
spond to the linear propagation (and mirror) phase shifts, 
modulo 27r. The A's, also called cavity detuning param- 
eters, are equal to {v — Uc)t, e.g. to the wave frequency 
mismatches from the nearest cold cavity frequency, scaled 
to the free spectral range l/r. 

The dynamics of the systems can be numerically stud- 
ied without any approximation by solving Eas. Hll3|l using 
a fourth-order Runge-Kutta solver with the appropriate 
initial conditions and making use of the boundary condi- 
tions Q (propagation model). Unstable fixed point can- 
not be found numerically so that approximate analytic 
solutions of the steady state equations must be worked 
out by expanding the amplitudes in Mac Laurin series 
of z, e.g. A,iz) = AjiO) + En=i[9i"^^i]oz"/n\, which 
allows to integrate Eas. ljll^ll . Such an expansion is justi- 
fied by the smallness of the scaled amplitudes {\Aj \ < 1) 
since giLi <C 1 and because z,z' < 1. The n-th order 
derivatives can be evaluated in terms of the field prod- 
ucts at z = using the generic equations O"®- To get 
the leading(fourth)-order coupling terms from this per- 
turbative approach, only the n = 1 terms in the field 
expansions need to be kept. After some algebra, the ap- 
proximate solutions of a-re 

Ap{t,Li + L2) = Ap{t,0) - {i/2)Ap{t,0)\As{t,0)\\{5) 

As{t,Li + L2) = A,{t,0) + {i/2)As{t,0)\Ap{t,0f 

- ^x*Al{t,0)[Alit,0)r, (6) 

A,(t,Li + L2) - iApit,0)A:it,0)+xA;{t,0)A^^it,m 

where the nonlinear coupling parameter is x = 
S'exp (— i^)(sin^/^). These solutions which assume a lin- 
ear z-variation of the fields depart from the numerical 
ones for decreasing pump resonance (SRO limit), but still 
account for pump depletion to first order. The cubic term 
in Eq.Q, present in conventional PRSRO model j^J, is 
due to the usual cascading (3a; — 2uj ~ iv) followed by 
the re-combination process (uj + 2uj ~ 3uj), while the last 
quartic term in Eq. (|SJ| , describing the two-step processes 
(3lu~2uj — u) followed by (lo^lo — 2lo\ leads to injection- 
locking. Note that this term is quadratic in the signal- 
and-idler resonant cases (see Eqs.(9) in Ref. 01) or in 
the pure IRO case without pump depletion (see Eq.(ll) 
of Ref. |5|). Stationary solutions to Eqs.Q)-© are ob- 
tained by requiring that v4p,s(t-f-T, 0) — ylp,s(t, 0) = Ap^s- 
Considering small enough detuning (Ap^s <C 27r), the ex- 
ponential phase factor in Eq.l^J is expanded as ~ l-l-iAj. 
The resulting steady state amplitude equations are 



(Kp-irpl\p)Ap = -(l/2)rp(l + iAp)Ap 1^1 + A, 

Ka iAe 



The above Mac-Laurin solutions (as compared with the 
mean-field approach based on amplitude expansion in the 
power of the k's 0) converge satisfactorily to the full 
propagation model as long as Vp > 0.8. Excellent conver- 
gence (to ±3% for the parameters of FiglSJ for instance) 
is found when both the pump and signal experience low 
roundtrip loss, not exceeding a few percent even for large 
pumping (/;„ = I Amp up to ~50 times the threshold for 
oscillation). For decreasing pump resonance and mod- 
erate pumping they still provide a qualitative account of 
the dynamical behavior of the system, although resulting 
in higher values of the intensities Ij — \Aj\^ . The domain 
of validity of the Mc Laurin model will be shown to de- 
pend on the value of Vp. In the true SRO Hmit (rp 0), 
the full propagation model remains always valid. 

Besides the trivial (non-lasing) solutions Ag^i — 0, 
Ap — A'lji I (^hip ^Tp Ap), Eqs.(|HJ admit non-zero inten- 
sity states Ip^s = l^p.sP with well-defined phases. It is 
convenient to introduce the scaled intensities Ip = Ip/I^ ^ 
Is = Is/Is, h = h/ill^Is), with = 2ks, Is = 2Kp/rp- 
input pump intensity lin = Iin/{2KpKs) ; and normal- 
ized cavity detuning Ap^s = Ap^s/ i^p,s- Introducing 
Cs ~ 4|x|^ (Kp/vp), and taking the moduhi of 0, the 
scaled signal intensity J, is then the solution of 



[F{l-rsKsAl) 



rs{l + nlAl)I,^Y + {AsFf 
^2CsrUl + KlAl)'lJ, 



2 

in ' 



where the symbol F stands for 



F=il + Ff + Al {KpF-rpf 



(9) 



(10) 



Tsil + iAs 



-As 



[MmAsApl^ ^ix*A'p[A:Y 



(8) 



The intracavity pump is given by FIp ~ /j„ and the idler 
intensity by It = /p/s[l + IslxlF + (Ip - l)/4]- 

Phase relationships, demonstrating phase-locking of 
the subharmonic waves to the pump laser, can be de- 
rived from Eqs.® by writing Aj = ctj exp (i(^j) where 
aj are the amplitude modulii. Defining <pi„ as the arbi- 
trary phase of the pump laser and (po = + 2ipp — Stps, 
one obtains 

tan((^m - ipp) = -[Ap(rp - Kpls)]/[l + F], (11) 
cot(pD ^ As/[1 - rsKsAl - rs{l + KlAj)Ip]. (12) 

When solved for (ps , these relations yield 3 possible values 
= fo + 2kTT/3 {ipo being a constant and k = 0,±1), 
while for \x\ = only the sum phase ips + fi = ipp + tt/2 
is deterministic as predicted for conventional oscillators, 
due to the phase diffusion noise stemming from the spon- 
taneous parametric fluorescence 0. Furthermore, when 
Ixl = (Cs = 0), Eq.Q imphes that the signal resonates 
necessarily with zero detuning (A^ = 0), and one re- 
trieves the result that the intracavity pump is clamped to 
the constant value Tp = 1 for any li^ O- From Ea. (|10|l . 
F is then the solution of a quadratic equation which ad- 
mits a single positive solution (supercritical bifurcation) 
if and only if /i„(l-t-KpAp) — Ap > 0. This condition is al- 
ways satisfied for > Ith , where the input (intracavity) 
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FIG. 1: Bifurcation diagram of signal intensity versus pump 
parameter, computed from the Mc Laurin solutions, for — 
0.3, Ap = 0, Ks = 0.005, X = S ^ 0.2. Curves (l)-(3) are for 
As = 0; 0.4; 0.5. The inset plots show the intracavity pump 
stable (solid line) and unstable (dashed) fixed points I^. Note 
that diverges at /i„ > 3.8; 4.4; 4.7 for curves 1;2;3 from 
either the LSA analysis or the time mapping of Eqs.l|2J-l|7J, 
setting the validity range of the Mc Laurin approximation. 



pump threshold expresses as Ith = l + '^pAp. Considering 
now the case |x| ^ 0, the signal wave is a priori no longer 
constrained to oscillate with zero detuning and its inten- 
sity is the solution of a quartic equation, X]ri=o ~ 0; 
obtained by expanding Eq.(|5J) using Ea. (|lUf) . The numer- 
ical resolution of this equation for a wide range of signal 
detuning or driving pump intensity always yield two real 
positive roots I^, defining two branches of solutions, for 
pump intensities > Iq (see below for definition of 
Jo). The stability of these two fixed points, each asso- 
ciated with the 3 possible phase states, was investigated 
using a linear stability analysis (LSA) that leads to 
a quartic characteristic equation, X]n=o ^n^"' — 0. The 
LSA results were double-checked by a direct time map- 
ping of Eqs.ljSl)-® using the boundary conditions 10} . 
From the LSA of the trivial state, the threshold for os- 
cillation (not necessarily on a stable fixed point) remains 
the same as for conventional devices. Fig^ shows the 
analytical bifurcation plot versus the pump parameter 
for Tp = 0.3 and 3 different values for Ag, as compared 
with the solutions computed with the propagation model 
(Fig|21). The bifurcation diagram versus /i„ displays a 
saddle-node region, with Iq being the input intensity at 
the saddle- node point where — 17- For zero or van- 
ishingly small As, one has Iq < Ith but as the detun- 
ing is increased, Iq > Ith- This behavior contrasts with 
the signal-and-idler resonant cases, for which a transition 
from sub-criticality to super-criticality is predicted when 
As,i 0, corresponding to the merging of the saddle- 
node point of coordinate {Iq,Is 0) with the threshold 
point {Ith ,Is =0) 3] . Note also the tiny range of sub- 
threshold states in curve (1), resulting from the weak 
self-injection regime. From the LSA of the Mc Laurin 
solutions, the I~ branches (dashed curves) are always 
unstable, due to a positive real eigenvalue, but only a 
portion of the upper branch (solid lines) extending from 
the saddle-node intensity Iq to some critical intensity Ic 
is found stable. The instability (dotted lines) beyond Ic 
(characterized by 3 real positive eigenvalues of the LSA 
equation) was confirmed by the time mapping of the Mc 
Laurin solutions which diverges at Ic, wherewhile the 



FIG. 2: Stable (upper) branches of signal intensity versus 
pump parameter, computed from the propagation model with 
the same parameters as for FigQ The hatched domain at the 
left of curve (2) gives the amplitude of the limit cycles below 
the saddle- node intensity Iq (non-stationary SPL states). The 
inset plots show the output {z' = 1) intracavity pump Ip. 



full propagation model converge to a fixed point for any 
input intensity value (Fig. |2J|. Hence the whole upper 
branch of the SPL-PRSRO is actually stable, the insta- 
bility predicted by the Mc Laurin model being merely an 
artefact of the approximation. Actually, as the pump en- 
hancement is decreased the validity range of the Mc Lau- 
rin model is restricted to pump parameters lying closer 
and closer to the threshold. In Fig. |2 obtained by back- 
ward adiabatic following of the stationary solutions, all 
curves end at their saddle-node point Iq since critical 
slowing down is observed as these points are approached. 
The limit cycles occurring for /i„ < Iq (hatched area) 
are merely due to the non-existence of stationary states. 
The upper branch is found stable even in the SRO limit 
(rp = 0), whatever the detunings, in contrast with the 
SPL-IRO mean-field analysis results Let us remind 
that in the SPL-DRO/TROs the upper branch was found 
to destabilize via a Hopf bifurcation ^] . The difference in 
dynamical behavior is due to the stronger self-injection 
regime in these latter devices. The inset frame in Fig. |21 
gives the output pump intensity (at z' = 1) correspond- 
ing to Ap s ~ (curve (1)), it can be seen that the 
pump is no longer clamped to unity as in conventional 
PRSROs meaning that the competing nonlinearity 
enhances the down-conversion efficiency. In Fig. |3| as 
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pump input I. 



FIG. 3: Bifurcation diagram of signal intensity versus pump 
parameter, computed from the analytical model for Vp = 0.95 
and same other parameters as in Fig0 excepted that curves 
(a)-(b) are for As = 0;0.1. The thin solid lines under-riding 
the (a)-(b) solid lines show the upper branch computed from 
the propagation model. The insets (A: Mc Laurin model, B: 
propagation model) are blow-ups of the saddle-node region 
of (b) case. In (B), critical slowing-down characteristic of 
saddle-node points occurs. 
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coupling strength S 



FIG. 4: Critical values of the signal detuning parameter 
delimiting the boundary of stable (underneath the curves) 
and oscillatory SPL states (above) in the (As, 5") plane at 
fixed pump input /detuning, computed from the propaga- 
tion model. The thick solid lines are for rp — 0.3 and the 
thin lines for rp — 0.95. The curve labels (l)-(3) stand for 
5 = 0; 7r/2, 37r/4. 




pump input I. 



FIG. 5: Domain of stability (dark region) of phase-locked 
states in the (5, /in) plane (propagation model), with Vp — 
0.95, Ks = 0.005, Ap = 0: (a),(c) for S = 0.2 and (b),(d) for 
S = 0.75. Upper (lower) frames are for Aa = (Aa = 0.05). 
The solid lines give the boundary of existence of stationary 
states from the Mc Laurin model. The plots are symmetric 
for negative phase mismatch ^. 

the pump finesse increases (np = 0.05), the Mc-Laurin 
(thick solid hnes) and the propagation model (thin solid 
hnes) converge excellently for any input intensity. Both 
models predict then a whole upper branch stability, for 
a wide range of input intensity (up to tested Im — 50 
at least). The Mc Laurin model converges to the prop- 
agation model because when Kp — > the assumption of 
a linear z-dependence of the resonating fields inside the 
medium is fully justified (uniform field limit). Notice 
that the threshold value is then more sensitive to the sig- 
nal detuning than in Fig. |21 As a consequence of the 
reduced sub-threshold state range due to the weak SPL 
regime one expect that the self-locking detuning range. 



defined as the maximum allowed for a given pump 
input intensity, would be smaller than in signal-and-idler 
resonant set-ups. Indeed large pump enhancement is 
paid back with a shrinking self- locking range. Fig0| dis- 
plays the critical Ac(S') detuning for both Vp — 0.3 and 
Tp = 0.95, when Ap = and when the device is pumped 
4 times above threshold. Surprisingly, for the same non- 
linear coupling strength S, the high pump resonance case 
leads to less than a cavity linewidth SPL range, while for 
the SRO limiting case the SPL range is 10-fold wider, 
as with SPL-DRO/TROs. Hence a large SPL range 
is not necessarily associated with the double-resonance 
condition that involves a strong mutual-injection process 
between the subharmonics. The small (< 1-MHz, e.g. 
less than a cavity linewidth) self-locking range reported 
by Boiler et al jlj] for a SPL-PRSRO is in agreement 
with this expectation. As a noticeable difference with 
the signal-and-idler resonant cases, the coefficients of the 
LSA equation depend on the phase of the coupling pa- 
rameter X, while X enters only as its modulus in the SPL- 
DRO stability analysis (see the Appendix in Ref.'3|). 
One hence expects that the dynamics of SPL-PRSROs 
will be sensitive to the SHG phase mismatch ^. FigEl 
shows the steady-state signal intensity contour plots in 
the parameter space (/i„,^), computed with the propa- 
gation model for /i„ 4, rp = 0.95, Ap = and two 
values S — 0.2 (a,c) and S = 0.75 (b,d) of the couphng 
parameter. For the stronger coupling S — 0.75 (panels 
b-d), small amplitude limit cycles arise around ^ = ±7r, 
even for (5s = (panel b). To avoid such non-stationary 
states, it is important to tailor accurately the simultane- 
ous phase-matching of both competing processes. 

In conclusion, it is found that the whole upper branch 
of the SPL- (PR) SRO is stable and that the self-locking 
range shrinks with increasing pump enhancement. In the 
SRO limit, the approximated model may fail to describe 
correctly the dynamical behavior over a large pump input 
range. Of more concern, in contrast with signal-and-idler 
resonant devices, the domain of existence of stationary 
phase-locked states is found sensitive to the value of the 
residual phase mismatch of the competing SHG nonlin- 
earity. The author is indebted to one of the referees for 
his personal involvement in the improvement of this re- 
port. This work has benefited from a partial support 
from an European Union INCO-Copernicus grant (Con- 
tract No. ERBIC15CT980814). 
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